Profiling diverse sequence tandem repeats in colorectal cancer reveals co-occurrence of microsatellite and chromosomal instability involving Chromosome 8.

We developed a sensitive sequencing approach that simultaneously profiles microsatellite instability, chromosomal instability, and subclonal structure in cancer. We assessed diverse repeat motifs across 225 microsatellites on colorectal carcinomas. Our study identified elevated alterations at both selected tetranucleotide and conventional mononucleotide repeats. Many colorectal carcinomas had a mix of genomic instability states that are normally considered exclusive. An MSH3 mutation may have contributed to the mixed states. Increased copy number of chromosome arm 8q was most prevalent among tumors with microsatellite instability, including a case of translocation involving 8q. Subclonal analysis identified co-occurring driver mutations previously known to be exclusive.

A genome-wide approach for detecting novel insertion-deletion variants of mid-range size.

We present SWAN, a statistical framework for robust detection of genomic structural variants in next-generation sequencing data and an analysis of mid-range size insertion and deletions (<10 Kb) for whole genome analysis and DNA mixtures. To identify these mid-range size events, SWAN collectively uses information from read-pair, read-depth and one end mapped reads through statistical likelihoods based on Poisson field models. SWAN also uses soft-clip/split read remapping to supplement the likelihood analysis and determine variant boundaries. The accuracy of SWAN is demonstrated by in silico spike-ins and by identification of known variants in the NA12878 genome. We used SWAN to identify a series of novel set of mid-range insertion/deletion detection that were confirmed by targeted deep re-sequencing. An R package implementation of SWAN is open source and freely available.

Haplotyping germline and cancer genomes with high-throughput linked-read sequencing.

Haplotyping of human chromosomes is a prerequisite for cataloguing the full repertoire of genetic variation. We present a microfluidics-based, linked-read sequencing technology that can phase and haplotype germline and cancer genomes using nanograms of input DNA. This high-throughput platform prepares barcoded libraries for short-read sequencing and computationally reconstructs long-range haplotype and structural variant information. We generate haplotype blocks in a nuclear trio that are concordant with expected inheritance patterns and phase a set of structural variants. We also resolve the structure of the EML4-ALK gene fusion in the NCI-H2228 cancer cell line using phased exome sequencing. Finally, we assign genetic aberrations to specific megabase-scale haplotypes generated from whole-genome sequencing of a primary colorectal adenocarcinoma. This approach resolves haplotype information using up to 100 times less genomic DNA than some methods and enables the accurate detection of structural variants.

Sensing of latent EBV infection through exosomal transfer of 5'pppRNA.

Complex interactions between DNA herpesviruses and host factors determine the establishment of a life-long asymptomatic latent infection. The lymphotropic Epstein-Barr virus (EBV) seems to avoid recognition by innate sensors despite massive transcription of immunostimulatory small RNAs (EBV-EBERs). Here we demonstrate that in latently infected B cells, EBER1 transcripts interact with the lupus antigen (La) ribonucleoprotein, avoiding cytoplasmic RNA sensors. However, in coculture experiments we observed that latent-infected cells trigger antiviral immunity in dendritic cells (DCs) through selective release and transfer of RNA via exosomes. In ex vivo tonsillar cultures, we observed that EBER1-loaded exosomes are preferentially captured and internalized by human plasmacytoid DCs (pDCs) that express the TIM1 phosphatidylserine receptor, a known viral- and exosomal target. Using an EBER-deficient EBV strain, enzymatic removal of 5'ppp, in vitro transcripts, and coculture experiments, we established that 5'pppEBER1 transfer via exosomes drives antiviral immunity in nonpermissive DCs. Lupus erythematosus patients suffer from elevated EBV load and activated antiviral immunity, in particular in skin lesions that are infiltrated with pDCs. We detected high levels of EBER1 RNA in such skin lesions, as well as EBV-microRNAs, but no intact EBV-DNA, linking non-cell-autonomous EBER1 presence with skin inflammation in predisposed individuals. Collectively, our studies indicate that virus-modified exosomes have a physiological role in the host-pathogen stand-off and may promote inflammatory disease.

Emergence of Hemagglutinin Mutations During the Course of Influenza Infection.

Influenza remains a significant cause of disease mortality. The ongoing threat of influenza infection is partly attributable to the emergence of new mutations in the influenza genome. Among the influenza viral gene products, the hemagglutinin (HA) glycoprotein plays a critical role in influenza pathogenesis, is the target for vaccines and accumulates new mutations that may alter the efficacy of immunization. To study the emergence of HA mutations during the course of infection, we employed a deep-targeted sequencing method. We used samples from 17 patients with active H1N1 or H3N2 influenza infections. These patients were not treated with antivirals. In addition, we had samples from five patients who were analyzed longitudinally. Thus, we determined the quantitative changes in the fractional representation of HA mutations during the course of infection. Across individuals in the study, a series of novel HA mutations directly altered the HA coding sequence were identified. Serial viral sampling revealed HA mutations that either were stable, expanded or were reduced in representation during the course of the infection. Overall, we demonstrated the emergence of unique mutations specific to an infected individual and temporal genetic variation during infection.

Metastatic tumor evolution and organoid modeling implicate TGFBR2 as a cancer driver in diffuse gastric cancer.

BACKGROUND: Gastric cancer is the second-leading cause of global cancer deaths, with metastatic disease representing the primary cause of mortality. To identify candidate drivers involved in oncogenesis and tumor evolution, we conduct an extensive genome sequencing analysis of metastatic progression in a diffuse gastric cancer. This involves a comparison between a primary tumor from a hereditary diffuse gastric cancer syndrome proband and its recurrence as an ovarian metastasis. RESULTS: Both the primary tumor and ovarian metastasis have common biallelic loss-of-function of both the CDH1 and TP53 tumor suppressors, indicating a common genetic origin. While the primary tumor exhibits amplification of the Fibroblast growth factor receptor 2 (FGFR2) gene, the metastasis notably lacks FGFR2 amplification but rather possesses unique biallelic alterations of Transforming growth factor-beta receptor 2 (TGFBR2), indicating the divergent in vivo evolution of a TGFBR2-mutant metastatic clonal population in this patient. As TGFBR2 mutations have not previously been functionally validated in gastric cancer, we modeled the metastatic potential of TGFBR2 loss in a murine three-dimensional primary gastric organoid culture. The Tgfbr2 shRNA knockdown within Cdh1-/-; Tp53-/- organoids generates invasion in vitro and robust metastatic tumorigenicity in vivo, confirming Tgfbr2 metastasis suppressor activity. CONCLUSIONS: We document the metastatic differentiation and genetic heterogeneity of diffuse gastric cancer and reveal the potential metastatic role of TGFBR2 loss-of-function. In support of this study, we apply a murine primary organoid culture method capable of recapitulating in vivo metastatic gastric cancer. Overall, we describe an integrated approach to identify and functionally validate putative cancer drivers involved in metastasis.

A programmable method for massively parallel targeted sequencing.

We have developed a targeted resequencing approach referred to as Oligonucleotide-Selective Sequencing. In this study, we report a series of significant improvements and novel applications of this method whereby the surface of a sequencing flow cell is modified in situ to capture specific genomic regions of interest from a sample and then sequenced. These improvements include a fully automated targeted sequencing platform through the use of a standard Illumina cBot fluidics station. Targeting optimization increased the yield of total on-target sequencing data 2-fold compared to the previous iteration, while simultaneously increasing the percentage of reads that could be mapped to the human genome. The described assays cover up to 1421 genes with a total coverage of 5.5 Megabases (Mb). We demonstrate a 10-fold abundance uniformity of greater than 90% in 1 log distance from the median and a targeting rate of up to 95%. We also sequenced continuous genomic loci up to 1.5 Mb while simultaneously genotyping SNPs and genes. Variants with low minor allele fraction were sensitively detected at levels of 5%. Finally, we determined the exact breakpoint sequence of cancer rearrangements. Overall, this approach has high performance for selective sequencing of genome targets, configuration flexibility and variant calling accuracy.

RVD: a command-line program for ultrasensitive rare single nucleotide variant detection using targeted next-generation DNA resequencing.

BACKGROUND: Rare single nucleotide variants play an important role in genetic diversity and heterogeneity of specific human disease. For example, an individual clinical sample can harbor rare mutations at minor frequencies. Genetic diversity within an individual clinical sample is oftentimes reflected in rare mutations. Therefore, detecting rare variants prior to treatment may prove to be a useful predictor for therapeutic response. Current rare variant detection algorithms using next generation DNA sequencing are limited by inherent sequencing error rate and platform availability. FINDINGS: Here we describe an optimized implementation of a rare variant detection algorithm called RVD for use in targeted gene resequencing. RVD is available both as a command-line program and for use in MATLAB and estimates context-specific error using a beta-binomial model to call variants with minor allele frequency (MAF) as low as 0.1%. We show that RVD accepts standard BAM formatted sequence files. We tested RVD analysis on multiple Illumina sequencing platforms, among the most widely used DNA sequencing platforms. CONCLUSIONS: RVD meets a growing need for highly sensitive and specific tools for variant detection. To demonstrate the usefulness of RVD, we carried out a thorough analysis of the software's performance on synthetic and clinical virus samples sequenced on both an Illumina GAIIx and a MiSeq. We expect RVD can improve understanding the genetics and treatment of common viral diseases including influenza. RVD is available at the following URL:http://dna-discovery.stanford.edu/software/rvd/.

Hierarchy in somatic mutations arising during genomic evolution and progression of follicular lymphoma.

Follicular lymphoma (FL) is currently incurable using conventional chemotherapy or immunotherapy regimes, compelling new strategies. Advances in high-throughput sequencing technologies that can reveal oncogenic pathways have stimulated interest in tailoring therapies toward actionable somatic mutations. However, for mutation-directed therapies to be most effective, the mutations must be uniformly present in evolved tumor cells as well as in the self-renewing tumor-cell precursors. Here, we show striking intratumoral clonal diversity within FL tumors in the representation of mutations in the majority of genes as revealed by whole exome sequencing of subpopulations. This diversity captures a clonal hierarchy, resolved using immunoglobulin somatic mutations and IGH-BCL2 translocations as a frame of reference and by comparing diagnosis and relapse tumor pairs, allowing us to distinguish early versus late genetic eventsduring lymphomagenesis. We provide evidence that IGH-BCL2 translocations and CREBBP mutations are early events, whereas MLL2 and TNFRSF14 mutations probably represent late events during disease evolution. These observations provide insight into which of the genetic lesions represent suitable candidates for targeted therapies.

A novel persistence associated EBV miRNA expression profile is disrupted in neoplasia.

We have performed the first extensive profiling of Epstein-Barr virus (EBV) miRNAs on in vivo derived normal and neoplastic infected tissues. We describe a unique pattern of viral miRNA expression by normal infected cells in vivo expressing restricted viral latency programs (germinal center: Latency II and memory B: Latency I/0). This includes the complete absence of 15 of the 34 miRNAs profiled. These consist of 12 BART miRNAs (including approximately half of Cluster 2) and 3 of the 4 BHRF1 miRNAs. All but 2 of these absent miRNAs become expressed during EBV driven growth (Latency III). Furthermore, EBV driven growth is accompanied by a 5-10 fold down regulation in the level of the BART miRNAs expressed in germinal center and memory B cells. Therefore, Latency III also expresses a unique pattern of viral miRNAs. We refer to the miRNAs that are specifically expressed in EBV driven growth as the Latency III associated miRNAs. In EBV associated tumors that employ Latency I or II (Burkitt's lymphoma, Hodgkin's disease, nasopharyngeal carcinoma and gastric carcinoma), the Latency III associated BART but not BHRF1 miRNAs are up regulated. Thus BART miRNA expression is deregulated in the EBV associated tumors. This is the first demonstration that Latency III specific genes (the Latency III associated BARTs) can be expressed in these tumors. The EBV associated tumors demonstrate very similar patterns of miRNA expression yet were readily distinguished when the expression data were analyzed either by heat-map/clustering or principal component analysis. Systematic analysis revealed that the information distinguishing the tumor types was redundant and distributed across all the miRNAs. This resembles "secret sharing" algorithms where information can be distributed among a large number of recipients in such a way that any combination of a small number of recipients is able to understand the message. Biologically, this may be a consequence of functional redundancy between the miRNAs.

LMP1 association with CD63 in endosomes and secretion via exosomes limits constitutive NF-κB activation.

The ubiquitous Epstein Barr virus (EBV) exploits human B-cell development to establish a persistent infection in ∼90% of the world population. Constitutive activation of NF-κB by the viral oncogene latent membrane protein 1 (LMP1) has an important role in persistence, but is a risk factor for EBV-associated lymphomas. Here, we demonstrate that endogenous LMP1 escapes degradation upon accumulation within intraluminal vesicles of multivesicular endosomes and secretion via exosomes. LMP1 associates and traffics with the intracellular tetraspanin CD63 into vesicles that lack MHC II and sustain low cholesterol levels, even in 'cholesterol-trapping' conditions. The lipid-raft anchoring sequence FWLY, nor ubiquitylation of the N-terminus, controls LMP1 sorting into exosomes. Rather, C-terminal modifications that retain LMP1 in Golgi compartments preclude assembly within CD63-enriched domains and/or exosomal discharge leading to NF-κB overstimulation. Interference through shRNAs further proved the antagonizing role of CD63 in LMP1-mediated signalling. Thus, LMP1 exploits CD63-enriched microdomains to restrain downstream NF-κB activation by promoting trafficking in the endosomal-exosomal pathway. CD63 is thus a critical mediator of LMP1 function in- and outside-infected (tumour) cells.

Exosomes: Fit to deliver small RNA.

Exosomes are specialized membranous nano-sized vesicles derived from endocytic compartments that are released by many cell types. Microvesicles are distinctive from exosomes in that they are produced by shedding of the plasmamembrane and usually larger in size (>1 µm). Exosome biogenesis involves the tightly controlled process of inward budding from the limiting membrane of multivesicular bodies (MVBs). This results in numerous intraluminal vesicles in the lumen of MVBs that contain distinct protein repertoires. It has been suggested that microvesicles shed by certain tumor cells hold functional messenger RNA (mRNA) that may promote tumor progression. We discovered that purified exosomes contain functional microRNAs (miRNAs) and small RNA, but detected little mRNA. Although a clear and decisive distinction between microvesicles and exosomes cannot be made and different subsets of exosomes exist, we speculate that exosomes are specialized in carrying small RNA including the class 22-25 nucleotide regulatory miRNAs. To demonstrate this we developed a co-culture system and found that exosomes are continuously secreted and transferred from Epstein Barr virus (EBV)-infected cells to uninfected neighboring cells. Throughout exosome transfer, the exogenous EBV-encoded miRNAs were delivered to subcellular sites of miRNA-mediated gene repression. Additionally, we found evidence that mature miRNAs are transferred between circulating cells in humans, since we detected EBV-miRNAs in non-infected cells in the peripheral blood of patients that include monocytes and T cells. In this addendum we discuss these findings in the context of recently published papers that advanced our current knowledge of exosome physiology, (mi)RNA function and intercellular RNA transfer. Based on this information we propose that an intercellular (miRNA-based) mode of signal transmission may be well suited in controlling space-confined processes such as the initiation of immune responses in the secondary (peripheral) lymphoid tissues or in a tumor microenvironment. Deciphering the molecular mechanism(s) that control small RNA loading into exosomes and transfer to recipient cells in vitro will provide new evidence for the physiological relevance of vesicle-mediated intercellular communication in vivo.

Functional delivery of viral miRNAs via exosomes.

Noncoding regulatory microRNAs (miRNAs) of cellular and viral origin control gene expression by repressing the translation of mRNAs into protein. Interestingly, miRNAs are secreted actively through small vesicles called "exosomes" that protect them from degradation by RNases, suggesting that these miRNAs may function outside the cell in which they were produced. Here we demonstrate that miRNAs secreted by EBV-infected cells are transferred to and act in uninfected recipient cells. Using a quantitative RT-PCR approach, we demonstrate that mature EBV-encoded miRNAs are secreted by EBV-infected B cells through exosomes. These EBV-miRNAs are functional because internalization of exosomes by MoDC results in a dose-dependent, miRNA-mediated repression of confirmed EBV target genes, including CXCL11/ITAC, an immunoregulatory gene down-regulated in primary EBV-associated lymphomas. We demonstrate that throughout coculture of EBV-infected B cells EBV-miRNAs accumulate in noninfected neighboring MoDC and show that this accumulation is mediated by transfer of exosomes. Thus, the exogenous EBV-miRNAs transferred through exosomes are delivered to subcellular sites of gene repression in recipient cells. Finally, we show in peripheral blood mononuclear cells from patients with increased EBV load that, although EBV DNA is restricted to the circulating B-cell population, EBV BART miRNAs are present in both B-cell and non-B-cell fractions, suggestive of miRNA transfer. Taken together our findings are consistent with miRNA-mediated gene silencing as a potential mechanism of intercellular communication between cells of the immune system that may be exploited by the persistent human gamma-herpesvirus EBV.

Human and mouse oligonucleotide-based array CGH.

Array-based comparative genomic hybridization is a high resolution method for measuring chromosomal copy number changes. Here we present a validated protocol using in-house spotted oligonucleotide libraries for array comparative genomic hybridization (CGH). This oligo array CGH platform yields reproducible results and is capable of detecting single copy gains, multi-copy amplifications as well as homozygous and heterozygous deletions as small as 100 kb with high resolution. A human oligonucleotide library was printed on amine binding slides. Arrays were hybridized using a hybstation and analysed using BlueFuse feature extraction software, with >95% of spots passing quality control. The protocol allows as little as 300 ng of input DNA and a 90% reduction of Cot-1 DNA without compromising quality. High quality results have also been obtained with DNA from archival tissue. Finally, in addition to human oligo arrays, we have applied the protocol successfully to mouse oligo arrays. We believe that this oligo-based platform using 'off-the-shelf' oligo libraries provides an easy accessible alternative to BAC arrays for CGH, which is cost-effective, available at high resolution and easily implemented for any sequenced organism without compromising the quality of the results.